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Abstract 



We consider the energy flow into gaps between hard jets. It was previously 
believed that the accuracy of resummed predictions for such observables 
can be improved by employing the kt clustering procedure to define the gap 
energy in terms of a sum of energies of soft jets (rather than individual 
hadrons) in the gap. This significantly reduces the sensitivity to correlated 
soft large-angle radiation (non-global leading logs), numerically calculable 
only in the large A^^^ limit. While this is the case, as we demonstrate here, 
the use of kt clustering spoils the straightforward single-gluon Sudakov ex- 
ponentiation that multiplies the non-global resummation. We carry out an 
0{af) calculation of the leading single- logarithmic terms and identify the 
piece that is omitted by straightforward exponentiation. We compare our 
results with the full ^(a^) result from the program EVENT2 to confirm 
our conclusions. For e~^e~ 2 jets and DIS (1+1) jets one can numerically 
resum these additional contributions as we show, but for dijet photoproduc- 
tion and hadron-hadron processes further studies are needed. 



1 Introduction 



Energy flow into gaps between hard jets is a valuable source of information on many 
aspects of QCD. Since this radiation is typically soft, a perturbative calculation of the 
corresponding distribution contains large logarithms that need rcsummation. Compar- 
isons of resummed perturbative estimates with data then enable us to test and further 
our knowledge of soft QCD dynamics. 

Moreover the hadronisation corrections are large, making these spectra a useful 
testing ground for theoretical ideas about power corrections within say a dispersive 

model for the QCD coupling [1]. Additionally at hadron colliders the activity away 
from jets has traditionally been used to study the soft underlying event and to refine 
models thereof [2], which will be an important component of physics at the LHC. 

In the present letter we analyse the first aspect alone, that of perturbative rcsum- 
mation. This resummation is challenging since the leading single-logarithms (there are 
no coUinear enhancements due to the fact that we are away from the hard jets) are gen- 
erated both by direct angular-ordered emission into the gap by the hard parton (jet) 
system as well as multiple energy-ordered correlated emission by a complex ensemble 
of soft emissions outside the gap [3,4], in addition to the hard jets. This latter piece, 
which cannot yet be computed analytically and more worryingly beyond a large Nc 
approximation, is known as the non-global component of the result. The non-global 
term causes a much stronger suppression of the soft energy flows than that obtained by 
vetoing direct emission off hard partons, into the gap, the Sudakov or Bremsstrahlung 
component of the answer [4]. The flnal result, for the cross-section with gap energy less 
than Qci can be expressed as [3-5] 

E(g, Qn) = exp [-R{Q/Qn)] S (Q/Qn) • (1) 

In the above exp[— i?] is the Sudakov term obtained by exponentiating the single 
gluon contribution and accounting for gluon branching to reconstruct the scale of the 
running coupling, while S is the non-global part of the answer obtained by running 
a large N^. Monte-Carlo program that encodes soft evolution of a system of dipoles 
to single-log (SL) accuracy, for emission into the gap fl. Equivalently for such away- 
from-jet energy flows one has to solve numerically a non- linear evolution equation [5] 
obtained in the large Nc limit. 

Given that only large Nc approximations of the non-global component are calculable 
at present in conjunction with the fact that it dominates the full result at smaller 
Qn values, it is important to reduce the sensitivity of the measurement to this effect. 
One method of doing so is to study event-shape-energy-flow correlations [6, 7] where 
measuring an event shape V outside the gap at the same time as the energy flow Qn 
inside it leads to non-global logarithms in the ratio V/Qn. Thus choosing V ~ Qn the 
magnitude of the non-global effects is reduced. However this procedure is more complex 
to implement in the case of several hard partons (e.g jet production in hadron-hadron 
collisions) and is quite restrictive experimentally, amounting essentially to studying a 
different observable. 
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Another version of the measurement was suggested in Refs. [8,9], to reduce the 
impact of non-global logarithms. There it was shown that defining the energy as 
the sum of energies of soft (mini) jets inside the gap region significantly reduced the 
non-global component of the result. Defining the jets via a kt clustering procedure [10], 
which is also common practice experimentally [11,12], had the effect of pulling soft 
emissions out of the gap region by clustering with harder emissions outside. As was 
shown in the simple case of e^e~ 2 jets [8], it is still possible for emissions near the 
centre of the gap to escape clustering. However since these emissions are well separated 
in rapidity from their nearest neighbours outside the gap (to escape clustering) the 
magnitude of non-global effects is reduced. This is because the bulk of the non-global 
piece arises in the region when the emitted soft gluon does not have too large an opening 
angle relative to those involved in the emitting ensemble [4]. Forcing a relatively large 
opening- angle /rapidity separation between the softest gluon and the harder emitters, 
as is required to escape clustering, does reduce the size of the non global effects [8] . 

We would like to point out, however, that using kt clustering impacts the general 
form Eq. ^ and does not leave the primary emission Sudakov piece exp[— i?] unchanged 
as has been assumed till now [8,9]. In fact we find the exact calculation of R becomes 
non-trivial at higher orders since it depends at rfi^ order on the n gluon geometry and 
the use of the clustering algorithm. The departure of R from its naive one-gluon form 
starts with two gluons and the resulting piece does not have the properties of non-global 
logs , neither in the colour structure nor the dynamical properties. This conclusion is 
unfortunate especially in the case of dijet photoproduction [9] or hadron-hadron studies 
where the missing piece we compute at leading order here, will have a highly non-trivial 
colour structure since it reflects the colour flow of the primary emission piece, computed 
e.g in [9,13]. This would impact accurate theoretical studies of such observables even 
though the non-global part is reduced. 

The current letter is organised as follows. In the subsequent section we identify the 
problem with naive one-gluon exponentiation to obtain the supposed Sudakov part of 
the answer, with kt clustering. We compute the piece that will be missed by one-gluon 
exponentiation, at leading order 0{a1). We then provide the full answer up-to order 
a1 , for the leading logarithms, and show via comparisons with EVENT2 [14] that 
the extra piece we compute is needed to agree with fixed order estimates while the 
expansion to order a1 of the Sudakov term, as currently computed in the literature, is 
insufficient. We conclude by pointing out that while analytical control over the Sudakov 
term is lost, due to clustering, it is possible in the simpler cases e^e~ two jets and 
DIS (1+1) jets, to numerically compute the additional piece at all orders, with existing 
programs [8], as we show. For higher jet topologies such as those studied in Ref. [9] it 
may only be possible to numerically compute this term in the large limit, thereby 
reducing the accuracy of the theoretical results compared to current expectations. 
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2 Independent soft gluon emission 



We specialise to the process e~^e~ 2 jets purely for reasons of simplicity. We perform 
our calculation for emission into a rapidity slice of width Arj, centred at zero rapid- 
ity with respect to the jet (thrust) axis, although similar considerations apply to any 
geometry one may choose for the inter-jet region Q. We first consider the Sudakov 
exponent generated by assuming exponentiation of single-gluon emission into the gap 
n (c.f. Eq. (3.9) of Ref. [8]): 

Sn.p(t) = - ^cko = e-^^-*^^ (2) 

with t = fg^"^ ^as{kt) = 4;^ln Yr2A ' where we used the one-loop running coupling 
to obtain up-to SL accuracy and defined A = /?oQ^s(Q) In 

We now argue that the above form which exponentiates the single gluon (leading- 
order) term is not correct starting from two gluon level. ^ Consider two soft emissions ki 
and ^2 that are ordered in energy ui ^ uj2, with ui outside the gap and uj2 inside. The 
situation is reminiscent of the correlated or non-global configuration at leading order. 
However that part pertains to the CfCa correlated gluon emission piece of the matrix 
element that is free from collinear singularities. In the present case we shall examine 
instead the independent emission Cp part of the two gluon matrix element along with 
the corresponding virtual corrections. We have for the independent emission of two real 
gluons by a dipole ab [15]: 

M'{k,, k,) = ClW^,{k,)W^,{k,) = 4Cl ^"f^ (3) 

[aki)[hki) {ak2){bk2) 

where Wab{ki) represents the emission of ki off the hard dipole ab and similarly for ^2- 

Now we examine the region where the two real soft gluons ki and k2 are clustered 
by the jet algorithm. This happens when 

(r7i-r72)' + (0i-02)'<i?', (4) 

where rj and (f) denote as usual rapidity and azimuth of the partons measured with 
respect to the axis defined by the back-to-back partons a and b and R is the radius 
parameter, usually set equal to one in experiment [11,12]. 

Since ki is outside the gap , ^2 is clustered into ki and pulled out of the gap. 
This configuration then does not contribute to the gap energy distribution da/duj (see 
Eq. ((21)). However now let us take the situation where we have ki virtual and /c2 as a 
real emission. Then ^2 is not clustered away by the algorithm and this configuration 
contributes with weight [15]: 

M2(/Cl,,irtual, k2) = -ClWabiki)Wabik2), UJl > UJ2- (5) 



^By this we mean that even after accounting for non-global logarithms, the exponentiation of the 
single gluon result still does not capture completely the remaining leading logarithms. 
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We thus have complete cancellation between the purely real and real-virtual terms, 
Eqs. (jni) and © in the region where ^2 is not removed by clustering. However in the 
angular region mapped by Eq. (0} only the one-real-one-virtual term will contribute, 
since k2 is in the gap. 

This contribution can then be computed as below. The four-vectors involved are 



a = 0,0,-1) (6) 

b = 0,0,1) (7) 

ki = A;(i(coshr7i, 0, 1, sinh?7i) (8) 

/c2 = A;t2(coshr72, sin(0), cos(0), sinh?72) (9) 



(10) 

where we have exploited the freedom to set (pi = and Q is the e~^e~ centre-of-mass 
energy. We have also neglected the recoil of the hard partons a and b, against the soft 
emissions ki and k2, which is valid for our aim of extracting the leading logarithms. 

Then the leftover real-virtual two-loop contribution reads (we compute the coeffi- 
cient of (as/27r)^) 

^„ = IQCl / ^ / fx2 - / dr^J dA{R'-{v2-Vir -<!>'). 

J X2 J xi V Q J Jk^tn Jk2en 27r 

(11) 

The above equation requires some explanation. We have introduced the dimen- 
sionless scaled transverse momenta xi,2 = 2/cfi,t2/Q and restricted the region such that 
virtual emission ki is integrated outside the gap region while real emission k2 inside. 
We have also inserted a step function that ensures that we are integrating over the 
region of Eq. where the corresponding double real emissions would be clustered 
and ^2 would be pulled out of the gap. The additional step function involving X2, that 
constrains the gap energy, is the usual one that corresponds to computing the cross- 
section for events with energy in the gap greater than Qq. This, by unitarity, is trivially 
related to that for events with gap energy less than Q^. From the latter quantity the 
distribution is directly obtained by differentiation with respect to Qq. We have denoted 
this term (JP"^'^'^^ jg ^ second order in as piece that has the colour structure and 
matrix element for independent emission from the primary dipole ab. However it is not 
derived by expanding the standard Sudakov result to order and is a companion to 
the non-global correction term 5*2 of [8], but with different functional properties and 
colour structure. 

Performing the integration over in Eq. (|lip we get 

^primary _ ^^L^ / dumm{u, At]) V - u\ (12) 
TT Jo 

Choosing for instance values of the gap size Arj > R we get 

^pnmary ^ ^cWR\ (13) 
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with L = In ^ . Alternatively choosing a smaller gap Arj < R we get 



^primary ^ 16^2^2 [1 (^3 _ (^2 _ ^^7])^') + ^ ( — 
TT 3 2 y 2 




) 



) 



(14) 



It is clear that although this piece has the same colour structure as that for inde- 
pendent two gluon emission, it cannot arise from expanding the single-gluon generated 
Sudakov Eq. (j21). The expansion of the naive Sudakov would give a term independent 
of R and which goes as {Ar])^ at O {aD- 

In the following section we shall show that the expansion of the Sudakov Eq. Q 
needs to be supplemented with the results Eqs. (HSl) , (HH) as appropriate, as well as the 
correlated non-global Ci?C^a^L^, piece computed in [8], in order to agree with the full 
a^L^ result, for Y,{Q/Q^) generated by the program EVENT2. 



3 Full 0{al) result and comparison to EVENT2 



where we denoted = %^ 

An additional CfCaoIL"^ term is indeed the non-global term computed in Ref. [8]. 
We compute this piece numerically for different values of the parameters Arj and R 
and add it to the result from Eq. (fT3j) for comparison with the fixed order program 
EVENT2. For example, with R = 1 and At] = 1 one gets 5*2 = —1.24:9 CfCa where 
5*2 is the first coefficient of the non-global log contribution S' = 1 + J2n=2 ^nt"', with t 
defined as before. 

The comparison to EVENT2 for the distribution da/dL is shown in Fig. Qfor the 
Cpa^L term, with L = IuQq/Q for R = 1 and Ar] = 1, as well as Arj = 0.5. If all lead- 
ing (single) logarithms in the integrated quantity T,{Q,Qq) are correctly accounted for 
by the resummed result Eq. Q, we would expect the difference between the EVENT2 
results and the expansion to NLO of the resummation, to be a constant at small Qq 
corresponding to large (negative) L. As we see this is only the case when (7|"''™^^y jg 
included by adding it to the expansion of the Sudakov Eq. Q. We have considered 
different values of R and Ar; as mentioned, for example, in the caption for Fig. ^ 
The comparison for other colour channels CfCa and CfTrUj shows agreement with 
EVENT2 (see Fig. 12)) which reflects the fact that only the channel, corresponding 
to independent emission, is incorrectly described by Eq. 



First we expand the Sudakov exponent Eq. (j21) to 0{a1). The result is 




(15) 
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4 All orders and conclusions 



We conclude by pointing out that in the simple cases of e'^e~ two jets and DIS 
(1+1) jets, the additional terms we describe here can be accounted for numerically, at 
all orders. This is done by using the Monte-Carlo program for large Nc dipole evolution 
developed in [3] and implemented with the kt clustering in [8]. By demanding multiple 
emissions from the primary dipole alone and restoring the colour factor for independent 
emission by changing Ca/2 —^Cp one obtains the result for primary emissions only, in 
the full theory. 

Using this procedure we see that the primary emission result, with kt clustering, 
differs from the Sudakov result generated by single gluon exponentiation, which is valid 
in the unclustered case. In Fig. |21 we show the results we obtain for primary emission 
with clustering and the Sudakov (unclustered) result. The discrepancy grows with the 
single-log evolution variable t = In where A = (3o<^s{Q) In ^ ~ ^■'^^ 

we note an increase of around 30% on inclusion of the terms we describe, that start 
with (jp^^^^^y computed analytically here. 

We wish to clarify that the full result in the large A''^ approximation, including the 
effect we point out here, is readily obtained by the method described in Ref. [8] and in 
fact computed there. However its separation into primary and non-global components 
(and restoring the proper colour factors where possible) needs to be done with care, 
keeping in mind our findings. 

The procedure to generate the most accurate theoretical results for the e'^e~ — >• 
two jets and DIS (1+1) jets is the following. We take the results as generated by the 
code used for Ref. [8], for a given gap geometry. This is the full result in the large Nc 
approximation and we divide it by the result obtained using the same code for primary 
emissions alone (rather than dividing by the naive Sudakov result), which takes as the 
only source for emissions the original hard dipole, e.g. the outgoing qq pair in e~^e~ 
annihilation. The result of this division is the non-global piece in the large limit. We 
can then make use of the fact that one can easily compute the exact 0{al) non-global 
term with proper CfCa colour factor and parameterise the non-global Monte-Carlo 
result, as a function of t, in a form that retains the correct colour structure for the 
leading a'^ln^ Q/Qq non-global term (see e.g. Ref. [3]). This is the non-global result, 
with the large Nc approximation starting only from O (o;^ In^ Q/Qn) terms. The overall 
result is obtained by multiplying the resultant parameterised form with the primary 
emission result, as computed here with full colour factors. The large N^. approximation 
is then confined to the non-global term and starting from the next-to-leading such 
piece (5*3 in the notation of [8]). It is thus still an important finding that the non-global 
logarithms are reduced considerably by kt clustering as demonstrated in Ref. [8], since 
this potentially reduces the impact of unknown non-global terms beyond the large Nc 
approximation. However the correct procedure for identifying the primary and non- 
global pieces, pointed out here, must be accounted for while comparing to experimental 
data to enable accurate phenomenological studies. 

In the case of dijet photoproduction, studied e.g. in [9], and gaps between jets 
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in hadron-hadron processes, it is less straightforward to account for the missing inde- 
pendent emission terms we point out. They will have a complex colour structure and 
existing large numerical programs cannot be employed to generate the full answer 
beyond the large A^^ limit. This would mean that the accuracy of the resummed result 
is limited not just by the unknown beyond-large-A'^c non-global logs but similarly in 
the primary emission terms which are not reduced by the use of clustering. In these 
cases further studies are therefore required to account correctly for the missing primary 
emission terms before one can argue that use of the clustering method mitigates the un- 
certainty involved in the theoretical predictions, by reducing the non-global component 
significantly. This is currently work in progress [16]. 
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Figure 1: Comparison of the C'^a'^L term produced by EVENT2 and the analytical 
calculation (referd to as resummed since it is derived by expanding the naive Sudakov 
resummation to NLO) with and without Cf'^™'^'^'^. The figures are for R = 1 and 
At] = 1.0 (above) and At] = 0.5 (below). The agreement for the unclustered case is 
also shown for comparison. 
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Figure 2: Comparison of the CfCa (above) and CpUf a^L term (below) produced by 
EVENT2 and the expanded Sudakov result, supplemented with non-global logs for the 
CpCa term. The figures are for R = 1 and A?7 = 1.0 and and as we expect the difference 
between the exact and resummed result expanded to NLO is a constant at large L. 
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Figure 3: The results for the primary emission resummation with and without kt clus- 
tering for R = 1 , At] = 1, using an adaptation of the program used for Ref. [8]. As can 
be seen, the clustering affects the primary emission term and the effect for t = 0.25 is 
an increment of over 30 %. 
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